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Sphaleron rate at high temperature in 1+1 dimensions 

Jan Smit and Wai Hung Tang a * 
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Valckenierstraat 65, 1018 XE Amsterdam, the Netherlands 

We resolve the controversy in the high temperature behavior of the sphaleron rate in the abelian Higgs model 
in 1+1 dimensions. The T 2 behavior at intermediate lattice spacings is found to change into T 2 / 3 behavior in the 
continuum limit. The results are supported by analytic arguments that the classical approximation is good for 
this model. 



Sphaleron physics plays an important role in 
theories of baryogenesis. A simple but useful 
model is the abelian Higgs model in 1+1 dimen- 
sions, given by 
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Recall that in 1+1 dimensions y/\ and g have di- 
mension of mass. A toy model for the electroweak 
theory is obtained by coupling to fermions, such 
that in the quantum theory the fermion current is 
anomalous. Changes in fermion number are then 
proportional to changes in Chern-Simons number 

1 



C = — dxi Ai, 
27r Jo 

thereby mimicking the B+L violation in the Stan- 
dard Model. Here we assume space to be a circle 
of circumference L. 

The sphaleron rate (of fermion number viola- 
tion) r can be identified from the diffusion of 
Chern-Simons number 



r = i([c(t)-c(o)] 2 



t 



For relatively low temperatures this rate is ex- 
ponentially suppressed by the sphaleron barrier. 
Numerical simulations |l|-|| aggree with analytic 
results H in this regime. At high temperature the 
rate is not known analytically, but expected to be 
un-suppressed. For temperatures larger than any 
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mass scale one may naively expect on dimensional 
grounds 



cx T 



(1) 



Such behavior was indeed found by us numeri- 
cally and reported at LATTICE 94 J§. However, 
De Forcrand, Krasnitz and Potting jy gave a scal- 
ing argument that the behavior should instead be 
given by 



~ cx T 2 '\ 



(2) 



This behavior (g) can also be argued for as fol- 
lows. The simulations use the classical approxi- 
mation 

([C7(i)-C7(0)] 2 )« 

e- H «te>*V T [C(tp(t), Tr(i)) - C(tp, n)} 2 
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Here tp and 7r denote generic canonical variables 
and ip(t) and 7r(t) are solutions of the classi- 
cal Hamilton equations with initial conditions 
ip(0) = tp, 7r(0) = 7r. The effective hamiltonian 
H c g is approximated by its classical form. Rescal- 
ing tp — > tpy/T, n — > 7rVT produces T only in the 
combination AT and e 2 T. Since this combination 
has mass dimension three, the behavior (||) in the 
form T/L cx (AT) 2 / 3 appears natural. 

We analyzed the quality of the classical approx- 
imation in perturbation theory and found the fol- 
lowing favorable properties (in 1+1 dimensions!) 

§: 
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Figure 1. New data for F = T /(m^L) versus T' 2 
for the lattice spacings given in (||) . The straight 
lines are fits to the data. The lower three lines do 
not fit the data very well. 



• Correlation functions of the basic fields are 
finite. 

• The approximation becomes exact in the 
weak coupling/high temperature limit 

A = ir 2 \ll\ 2 , g 2 = t\, T = v 2 \li\T', 
v 2 ~* oc, £,T' fixed. (3) 

Here ft is the renormalized mass parameter 
in the MS scheme and v 2 = —fi 2 /\ is the 
classical ground state value of 2ip*(p. In the 
limit (§), [i 2 III 2 -> 1. Notice that T = 
AT/|/i| involves the combination AT. 

In our previous results H the coefficient of T 2 
appeared to vanish on extrapolation of the lat- 
tice distance to zero. At the time we interpreted 
this as an effect caused by using a too simple ef- 
fective hamiltonian (the classical one), but now 
the perturbative analysis tells us that using the 
classical hamiltonian is fine in the limit (|J). To 



settle the issue we carried out additional simula- 
tions at higher temperatures and smaller lattice 
spacings. 

Fig. 1 shows the dimensionless rate 



F 



mz.L 



plotted versus T' 2 , for several lattice spacings, 



a\n\ = 0.25, 0.23, 0.16, 0.11, 0.08. 



(4) 



There is clear T 2 behavior for the two larger spac- 
ings, but for the three smaller spacings this be- 
havior does not fit the data any more. A crossover 
appears to take place between a|/i| = 0.23 and 
0.16. In fact, T 2 / 3 behavior fits the data better 
at the three smaller spacings. Figs. 2-3 illustrate 
the behavior in more quantitative detail. Fitting 
the forms 



F 2 = 
F2/3 
led to 



CO 
Co 



c 2 T L 
c 2 / 3 



rp'2/3 



x|/d.o.f. = 
X 2 /3 /d.o.f. = 



1.4, 0.42, 6.5, 4.8, 3.45, 
28.9, 21.6, 2.5, 1.1, 1.2, 



for a\fi\ = 0.25, 0.23, 0.16, 0.11, 0.08, respectively. 
The volume was fixed at \n\L = 16, and £ = 4. 

Clearly, the T 2 / 3 form is favored for the three 
smaller lattice spacings. Extrapolating the re- 
sulting c 2 /3 to zero lattice spacing, assuming a 
quadratic dependence on a (cf. Fig. ||), gives the 
result c 2 /3 = 0.00452(86) for a — 0, which trans- 
lates into 

y = 0.0090(17)(AT) 2/3 , 5 2 /A = 4. 

Arguments for a T 2 / 3 behavior were given in 
ref. [0, where it was shown that the classical 
rate could be written in terms of a function 
g(x,y), as T/L = T 2 / 3 g(a 3 T,v 3 /T), using units 
in which A = 1, suppressing £ dependence. For 
high temperatures T 2 / 3 behavior then followed 
in the continuum limit provided g(0, 0) 7^ 0, 00 
existed. For the case v 2 = 0, numerical data sup- 
ported a nontrivial <?(0, 0). Our analytic results 
^ support the existence of a continuum limit 



3 



of g(a 3 T,v 3 /T), but notice how the combina- 
tion a 3 T implies non-commutativity of the limits 
a — > and T — > oo. Unfortunately, our attempts 
to extrapolate first to zero lattice spacing failed 
because the resulting errors got too large. Our 
present analysis seems to indicate a behavior like 
g(x, y) = g(Q, 0)+j3(xy) 2 / 3 +ax 4/3 +jy 4/3 +- ■ ■, or 
T/L = [g(0,0) + f3a 2 }T^ 3 + aa 4 T 2 + 1 T- 2 / 3 + -- ■. 
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Figure 2. Plot of c 2 /3 as a function of a 2 \n\ 
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Figure 3. Top: data showing T 2 / 3 behavior for 
a\n\ = 0.16. Plotted is F versus T' 2 / 3 and a fit to 
Co +C2/3T' 2 / 3 . Bottom: same data plotted versus 
T' 2 . The straight line is a fit of c + c 2 T' 2 . The 
curved line is the T' 2 ! 3 fit of the top figure. The 
data lack T 2 behavior but favor T' 2 / 3 behavior in 
the region a\fi\ < 0.16. 



